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We present an introduction to dynamical trapping horizons as quasi-local models for 
black hole horizons, from the perspective of an Initial Value Problem approach to the 
^ ' construction of generic black hole spacetimes. We focus on the geometric and structural 

properties of these horizons aiming, as a main application, at the numerical evolution 
■ and analysis of black hole spacetimes in astrophysical scenarios. In this setting, we dis- 

nJ ' cuss their dual role as an a priori ingredient in certain formulations of Einstein equations 

(N ■ and as an a posteriori tool for the diagnosis of dynamical black hole spacetimes. Com- 

plementary to the first-principles discussion of quasi-local horizon physics, we place an 
emphasis on the rigidity properties of these hypersurfaces and their role as privileged 
geometric probes into near-horizon strong-field spacetime dynamics. 
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^ 1. Black holes: global vs. (quasi-)local approaches 

1.1. Establishment's picture of the gravitational collapse 

Our discussion is framed in tire problem of gravitational collapse in General Relativ- 
ity. The current understanding is summarized in what one could call the establish- 
ment's picture of gravitational coUapsJI^, a heuristic chain of results and conjectures: 

(1) Singularity Theorems: if gravity is able to make all light rays locally coriverq e 
(namely, if trapped surfaces exist), then a spacetime singularity formfPEEET 

(2) (Weak) Cosmic Censorship (Conjecture): in order to preserve predictability, the 
formed singularity is not visible for a distant observer^. 

(3) Black hole spacetimes stability (Conjecture): General Relativity gravitational 
dynamics drives eventually the black hole spacetime to a stationary state. 

(4) Black Hole uniqueness theorem: the final state is a Kerr black hole spacetimJ^. 

Light bending is a manifestation of spacetime curvature and black holes constitute a 
dramatic extreme case of this. The standard picture of gravitational collapse above 
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suggests two (complementary) approaches to the characterization of black holes: 

a) Global approach: (weak) cosmic censorship suggests black holes as no-escape 
regions not extending to infinity. Its boundary defines the event horizon <§ . 

b) Quasi-local approach: singularity theorems suggest the characterization of a 
black hole as a spacetime trapped region where all light rays locally converge. 

The establishment 's picture of gravitational collapse depicts an intrinsically dynami- 
cal scenario. Hence, a systematic methodology to the study of dynamical spacctimes 
is needed. We adopt an Initial (Boundary) Value Problem approach, that offers a 
systematic avenue to the qualitative and quantitative aspects of generic spacetimes. 

1.2. The Black Hole region and the Event Horizon 

The traditionalSl approach to black holes involves global spacetime concepts, in 
particular a good control of the notion of infinity. Given a (strongly asymptotically 
predictable) spacetime Ai, the black hole region is defined a& ^ — J^(J^+), 
where J~(^J'^) is the causal past of future null infinity J^+. That is, SB is the 
spacetime region that cannot communicate with J^+. 

We are particularly interested in characterizing a notion of boundary surface of 
black holes. In this global context this is provided by the event horizon (f , defined as 
the boundary of that is = 9J^(J^+) OAi. Interesting geometric and physical 
properties of the event horizon are: i) S' is null hypersurface in Ai; ii) it satisfies an 
Area Theorert^^, so that the area of spatial sections SofS' does not decrease in the 
evolution; and, beyond that, iii) a set of black hole mechanics laws are fulfiUedlS. 

However, the global aspects of the event horizon also bring difficulties: a) it is 
a teleological concept, i.e. the knowledge of the full (future) spacetime is needed in 
order to locate S', and b) the black hole region and the event horizon can enter into 
flat spacetime regions. In sum, the notion of event horizon is a too global one: it 
does not fit properly into the adopted Initial Value Problem approach. 

1.3. The Trapped Region and the Trapping Boundary 

The global approach requires controlling structures that are not accessible during 
the evolution. In this context, the seminal notion of trapped surfac^plays a crucial 
role, capturing the idea that all light rays emitted from the surface locally converge. 
Through the singularity theorems and weak cosmic censorship, it offers a bench- 
mark for the existence of a black hole region: in strongly predictable spacetimes 
with proper energy conditions, trapped surfaces lie inside the black hole regiorPl. 
Moreover, their location does not involve a whole future spacetime development. 

1.3.1. Trapped and outer trapped surfaces. Apparent horizons 

Given a closed spatial surface S in the spacetime, we can consider the light emitted 
from it along outer and inner directions given, respectively, by null vectors £° and 
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/c°. Then, light locahy converges (in the future) S if the area of the emitted light- 
front spheres decreases in both directions (see though Ref . [TT|) . Denoting the area 
element S as dA = y/q(Px, the infinitesimal variations of the area along i'^ and 
define outgoing and ingoing expansions 9^^^ and 6'^'""' (see section [01 for details) 

A trapped surface is characterized by O'-^^O^''^ > 0. In the black hole context, in 
which the singularity occurs in the future, we refer to 5 as a future trapped surface 
(TS) if 6*^^) < 0, fl^'^-' < and as future marginally trapped surface (MTS) if one of the 
expansions, say 9^^\ vanishes: 6''^-' = 0,9^''^ < 0. If a notion of naturally expanding 
direction for the light rays exists (e.g. in isolated systems, the owier null direction i^^ 
pointing to infinity), a related notion of outer trapped surface is give iPby 6'(^) < 0. 
Marginally outer trapped surfaces (MOTS) are characterized by 9^^^ = 0. 

Before proceeding to a characterization of black holes in terms of trapped sur- 
faces, let us consider trapped surfaces from the perspective of a spatial slice of 
spacetime S. The trapped region in S, 7s C S, is the set of points p G Ti belong- 
ing to some (outer) trapped surface 5 C S. The Apparent Horizon (AH) is then 
the outermost boundar y of the trapped region Ts . A crucial result is the following 
characterizatiorPIEIlSlQf AHs: if the trapped region 7s in a slice S has the structure 
of a manifold with boundary, the Apparent Horizon is a MOTS, i.e. 6*^ = 0. 

Given a 3-1-1 foliation of spacetime {St}, let us consider the worldtube obtained 
by piling up the 2-dimensional AHs St C . Such an AH- worldtube does not need 
to be a smooth hypersurface (it is not even necessarily continuous, as discussed in 
section fS.l.ip . This is our first encounter with the notion of a spacetime worldtube 
foliated by MOTS. Though these worldtubes are slicing-dependent, their character- 
ization in terms of MOTSs makes them very useful from a operational perspective. 

1.3.2. The trapped region: definition and caveats 

From a spacetime perspective, no reference to a slice E must enter into the charac- 
terization of the trapped region. The spacetime trapped region T is defined as the 
set of points p E Ai belonging to some trapped surface S C Ai. Its boundary is 
referred^ to as the trapping boundary. These concepts offer, in principle, an intrin- 
sically quasi-local avenue to address the notion of black hole region and black hole 
horizon, with no reference to asymptotic quantities. 

In spite of their appealing features, there are also important caveats associated 
with the trapped region and the trapping boundary. In particular, we lack an op- 
erational characterization of the trapping boundary (see also the contribution by 
J.M.M. Scnovilla). A systematic attempt to address this issue is provided by the 
notion of trapping horizoi^^, namely smooth worldtubes of MOTS (see section [2?2| . 
as a model for the trapping boundary. Trapping horizons, that are non-unique, have 
led to important insights into the structure of the trapped region, though an oper- 
ational characterization of the trapping boundary is still missing. 



August 12, 2011 0:9 WSPC/INSTRUCTION FILE JaramilloShanghai 



4 J.L. JARAMILLO 

The difficulties are illustrated in the discussion of the relation between the trap- 
ping boundary and S . In strongly predictable spacctimes with appropriate energy 
conditions (see, though Ref iTSj) . the trapped region T is contained i n the b lack hole 
region SS. In attempts to refine this statement, support was suggest- 
ing that the trapping boundary actually coincides with the event horizon, though 
later wor liff^ showed that the trapped region not always extend up to (a'. The ques- 
tion is still open for (outer) trapped regions constructed on outer trapped surfaces, 
rathe r than , on TSs. Important insight into these issues has been gained in recent 
wor demonstrating truly global features of the trapped region T. In partic- 

ular: 

i) The trapping boundary cannot be foliated by MOTS. 

ii) Closed trapped surfaces can enter into the fiat region. This is an important 
issue in this approach to black holes, since it was a main criticism in ll.2l 

iii) Closed trapped surfaces are clairvoyant^ that is, they are aware of the geometry 
in non-causally connected spacetime regions. This non-local property challenges 
their applicability for an operational characterization of black holes. 



1.4. A pragmatic approach to quasi-local black hole horizons 

Trapping horizons offer a sound avenue towards the quasi-local understanding of 
black hole physics. They provide crucial insight in gravitational scenarios where a 
quasi- local notion of black hole horizon is essential, such as black hole thermodynam- 
ics beyond equilibrium, the characterization of physical parameters of strongly dy- 
namical astrophysical black holes (notably in numerical simulations) , semi-classical 
collapse, quantum gravity or mathematical relativity (cf. A. Nielsen's contribution). 
But, on the other hand, issues like their non-uniqueness or the clairvoyant properties 
of trapped surfaces pose fundamental questions that cannot be ignored. 

We do not aim here at addressing first-principles questions about the role of 
trapping horizons as a characterization of black hole horizons. We rather assume 
a pragmatic approach to the study of gravitational dynamics, which underlines the 
role of trapping horizons as hypersurfaces of remarkable geometric properties in 
black hole spacctimes. More specifically, our main interests are: 

i) The construction and diagnosis of black hole spacetimes in Initial (Boundary) 
Value Problem approach. 

ii) Identification of a geometric probe into near-horizon spacetime dynamics. 

Point ii) is particularly important in the study of gravity in the strong-field regime, 
where the lack of rigid structures (e.g. symmetries, a background spacetime...) is 
a generic and essential problem. Given our interests and the adopted pragmatic 
methodology, we look for a geometric object such that: a) represents a footprint of 
black holes, providing a probe into their geometry; b) is adapted, by construction, 
to an Initial-Boundary Value Problem approach; and c) although not-necessarily 
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unique, provides a geometric structure with some sort of rigidity property. As we 
shall see in the following, dynamical trapping horizons fulfill these requirements. 

1.5. General scheme 

In section[2]we introduce the basics of the geometry of closed surfaces in a Lorentzian 
manifold and motivate quasi-local horizons in stationary and dynamical regimes. 
Section [3] reviews their geometric properties and their special features as physical 
boundaries. Sections H] and [5] are devoted to applications in a 3+1 description of the 
spacetime. Section |4] shows the use of quasi- local horizons as inner boundary condi- 
tions for elliptic equations in General Relativity, whereas section [5] discusses some 
applications to the analysis of spacetimes, in particular their role in a correlation 
approach to spacetime dynamics. In section [6] a general overview is presented. 

2. Quasi-local horizons: Concepts and Definitions 
2.1. Geometry of spacelike closed 2-surfaces S 

2.1.1. Normal plane: outgoing and ingoing null vectors 

Let us consider a spacetime (A4,gab) with Levi-Civita connection Vq. Given a 
spacelike closed (compact without boundary) 2-surface S in Ai and a point p £ S, 
the tangent space splits as TpM = TpS®Tjj-S. We span the normal plane T^l^S either 
by (future-oriented) null vectors £° and fc" (defined by the intersection between T^S 
and the null cone at p) or by any pair of normal timelike vector n° and spacelike 
vector s°. Let us denote conventionally to be the outgoing null normal and fc" 
the ingoing one. We choose normalizations: 

= , fc^fca = , rfca = -1 , n'^Ua = -1 , .s".s, = 1 , n'^Sa = , (2) 

Directions i"^ and are uniquely determined, but a normalization-boost freedom 

e'" = /r , k'" ^ r^k" (3) 

n'" = cosh(fT)n'' + sinh(CT)s'' , s'" = sinh(cr)n° -|- cosh(cr)s° , 

remains for some arbitrary rescaling positive function / on S (where a = ln(/) and 
r = A(n'^ + s'')/\/2 and k" = \-^{n'' - s")/V2, for some function A on S). 

2.1.2. Intrinsic geometry of S 

The induced metric on S is given by 

Qab = gab + kah + iakb = gab + riaUb ~ SaSb , (4) 

so that is the projector onto S 

<fbi\^<fc , q\v'' ^v-{\fv- eTS) , q\w'' ^OiVw'^ eT^S) . (5) 
We denote the Levi-Civita connection associated with Qab as ^Da- The volume form 
on S will be denoted by \ — ^/qdx^ A dx^, i.e. ^ab = n'^s'^'^^cdab, though we will also 
employ the area measure notation dA = y^d^x. 
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2.1.3. Extrinsic geometry of S in {AA,g) 

We define the second fundamental tensor oi {S,qab) in {M,gab) (also, shape tensor 
or extrinsic curvature tensor) as 

ICU^q\q\^dq\ , (6) 

where c is an index in the normal plane T-^S, whereas a and b are indices in TS. 
Given a vector u° normal to S, wc can define the deformation tensor 0^^' as 

Qil^-<l\Q\^cV, . (7) 
Then, using expression Q , the second fundamental tensor can be expressed as 

icib = fc^eit' + tei'^ = n^e^i^ ~ s^e^:^ . (8) 

(v) 

We can express in terms of the variation of the intrinsic metric along 
Given a (tcnsorial) object Aq^ . q^''^ "^'" tangent to iS we denote by Sv the operator 
{Sv A) ai...a„''^ = qa^""^ ■■■qa^''''qdt^^ -qdj'"^ CyAc^,,,cJ^ , whcrc £„ dcuotc the 
Lie derivative along (some extension of) Then, it follows 

^vqab = ieil^ . (9) 



a) Shear and expansion associated with v"". Defining the expansion 9^"^ and shear 
tensor cr^^-* associated with the normal vector u° as 

q-by^v, ^ S^ln^q , ^i;' ^ e^^,^ - ^e^^Uab , (10) 



2 



we express the deformation tensor 0^^^ in terms of his trace and traceless parts 

1 
2 

b) Mean curvature vector H'^. Taking the trace of 0^^^' on S we define the mean 
curvature vectoM 



il' = -i;' + ^S^"'qab . (11) 



F= = (7''''/C;^b = 6'(^)fc^ + 6l('^)r . (12) 

The extrinsic curvature information of (S^qat) in {A4, gab) is completed by the 
normal fundamental forms associated with normal vectors u". In particulap^ 

^^^^ = ^k^q\Vdc , = -^^V.Vrffc, . (13) 

All these normal fundamental forms are related up to a sign and a total derivative 
on S. Using the normalizations {2]) wc gclE|: f^i"^ = -ni'\da'^ = ~n''a\fl''a^ = 
fii"^ — ^DaX. We choose to employ the 1-form fla ^ in the following. 

^Note the opposite sign convention with respect to the contribution by J.M.M. Senovilla. 
''When using ^ka = -e" one gets: f^i^' = -Q.^a^ - "^Da^- This wiU be relevant later, in Eq. l|39|l . 
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2.1.4. Transformation properties under null normal rescaling 

Under the rescaling ([2]) £° — > fc" — > f ^^k"- the introduced fields transform as 

T^lb ^ Kb ^ 

^'2^f^'2 ei^^^fei^^ -it^^Mt' (14) 

r!i^)^r!l^)+2i?,(ln/) 
Finally, given an axial Killing vector (/)° on 5, wc can write the angular momenturr@ 

J=±fni'}r'e. (15) 
The transformation rule of fii^' in ([T4| together with the divergence-free property 



of 0" (following from its Killing character) guarantee that the quantity J does not 
depend on the choice of null normals fc" (i.e. J does not change under a boost). 



2.2. Trapping Horizons 

2.2.1. Worldtubes of marginally trapped surfaces 

A trapping horizoi^^is (the closure of) a hypersurface H foliated by closed marginal 
(outer) trapped surfaces: = [J^^^St, with ^^^^j^j^ = 0. Trapping horizons are 
classified according to the signs of 0^'^-' and Skd'^^^ In particular, the sign of 9'^'^^ 
controls if the singularity occurs either in the future or in the past of S, whereas 
the sign of SkO^^'^ controls the (local) outer- or innermost character of H. Then, a 
trapping horizon is said to be: i) future (respectively, past) if 9'-'^'> < (respectively, 
Q{k) > Q)^ outer (respectively, inner) if there exist^ £° and k'^ such that 

46''^' < (respectively, 46''^' > 0). 



2.2.2. Future Outer Trapping Horizons 

In a black hole setting the singularity occurs in the future of sections St of "H, so 
that the related trapping horizon is of future type, 0^'^-' < 0. In addition, when 
considering displacements along (ingoing direction) we should move into the 
trapped region, i.e. (J^^^^' < 0, so that the trapping horizon should be outer. 

The resulting characterization of quasi-local black hole horizons as Future Outer 
Trapping Horizons (FOTHs) is further supported by the following analysis of the 
area evolution. Hawking's area theorem for event horizons (cf. section [T|2]) captures 



"^The quantity J coincides with the Komar angular momentum in case that tp'^ can be extended 
to an axial Killing in the neighbourhood of cS. 

"^The sign of Sk^^^"' is not invariant on the whole S under a rescaling ((2}. However, if there exists 
l'^ and fe" such that Sf^d^^") < on 5, then there does not exist any choice of and fc" such that 
Sf^d^^^ > on (S; see Ref. 1221 and also the marginally trapped surface stability condition in Ref. 1231 
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a fundamental feature of classical black holes. It is natural to wonder about a quasi- 
local version of it. Let us consider an evolution vector along the trapping horizon 
T-L, characterized as: i) h"- is tangent to T-L and orthogonal to St, and ii) /i° transports 
St onto St+st'- ^ht — 1. We can write ft." and a dual vector r" orthogonal to H as 

Then h°-ha = —r'^Ta = 2C, i.e. is spacelike for C > 0, null for C = and timelike 
for C < 0. The evolution of the area A — dA = \ along /i° is given by 

6nA = y = y - c^^f'^') = - / ce^'^'e . (17) 

Considering for simplicity the spherical symmetric case (C = const; see discussion 
of Eq. ([57)1 in 13. 2. 41 for the general case), the trapping horizon condition, 5hd^^^ = 0, 
writes SgO^^^ — C5kO^^^ = 0, so that C ~ |^|w- ^PPlyi^ig the Raychaudhuri equation 
for (J^e^^) [see later Eq. (pij) ]. together with the 6''^^ = condition, we find 

• 

Under the null energy and outer horizon conditions, it follows C > 0, so that 
the future condition guarantees the non-decrease of the area in (|17p . Therefore, 
FOTHs are null or spacelike hypersurfaces (C > 0), satisfying an area law result, 
and therefore providing appropriate models for quasi-local black hole horizons. 



2.3. Isolated and Dynamical Horizons 

The distinct geometric structure of null and spatial hypersurfaces suggests differ- 
ent strategies for the study of the stationary and dynamical regimes of quasi-local 
black holes, modeled as future outer trapping horizons. This has led to the parallel 
development of isolated horizon and the dynamical horizon framework j^^ l ^^ ^ ^^ l ^^ i 
In equilibrium. Isolated Horizons (IH) provide a hierarchy of geometric struc- 
tures constructed on a null hypersurface H that is foliated by closed (outer) 
marginally trapped surfaces. They characterize different levels of stationarity for 
a black hole horizon in an otherwise dynamical environment: 

i) Non-Expanding Horizons (NEH). They represent the minimal notion of equilib- 
rium by imposing the stationarity of the intrinsic geometry Qat- 

ii) Weakly Isolated Horizons (WIH). They are NEHs endowed with an additional 
structure needed for a Hamiltonian analysis of the horizon and its related (thermo- 
)dynamics. They impose no additional constraints on the geometry of the NEH. 

iii) Isolated Horizons (IH). These are WIHs whose extrinsic geometry is also invari- 
ant along the evolution. They provide the strongest stationarity notion on H. 

The non-stationary regime can be characterized by Dynamical Horizons (DH), 
namely spacelike hypersurfaces H foliated by closed future marginally trapped sur- 
faces, i.e. fl'^-* = and fl'*^-* < 0. Introduced in a 3+1 formulation, they provide 
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a complementary perspective to the dual-null foliation formulatioi^^ of trapping 
horizons, making them naturahy adapted for an Initial Value Problem perspective. 

2.4. IHs and DHs as stationary and dynamical sections of FOTHs 

A natural question when considering the transition from equilibrium to the dy- 
namical regime is whether a section St of a FOTH can be partially stationary and 
partially dynamical. Or, in other words, whether the element of area dA can be non- 
expanding (C = 0) in a part of St whereas it already expands (C > 0) in another 
part. Namely, can h°- be both null and spacelikc on a section St of a FOTH? 

The answer is in the negative. Transitions between non-expanding and dynam- 
ical parts of a FOTH must happen all at once. More precisely, assuming the null 
energy condition, a FOTH can be completely partitioned into non-expanding and 
dynamical sections. For a section St to be completely dynamical (C > 0) it s uffices 
that it has di0^^'> < somewhere on it. Otherwise /i" is null (C = 0) all over 5j 28 | 22 | 

In more physical terms, it suffices that some energy crosses the horizon some- 
where, and the whole horizon instantaneously grows as a whole. This non-local 
behaviour is a consequence of the elliptic nature of quasi-local horizons. As shown 
in section I3.2.3[ the function C determining the metric type of /i" satisfies an el- 
liptic equation [cf. Eq. ([57]) ]. Under the outer condition Skd'-^'' < one can apply 
a maximum principle to show that C is non- negative [generalization of Eq. psp ]. 
Moreover, it suffices that SgO'-^'' ^ somewhere, for having C > everywhere. 

3. Quasi-local horizons: properties from a 3+1 perspective 
3.1. Equilibrium regime 

3.1.1. Null hyper surf aces: characterization and basic elements 

A hypersurface % is null if and only if the induced metric is degenerated. Equiv- 
alently, if and only if there is a tangent null vector I"' orthogonal to all vectors 
tangent to U: l^Va = 0, Vu" e TU. 

Let us introduce some elements on the geometry of Choosing a null vector 
transverse to we can writcEl the degenerated metric as qah = 5ab + ka^b + ^akb- 
A projector onto H can also be constructed as: IlJ' — Sj' + £ak^ ~ qj' — h^i^ . As 
a part of the extrinsic curvature of "H, a rotation 1-form can be introduced^ on H 

it] 1 (£) 

as Wa = ■j^k'^^ a^c- This 1-form lives on i.e. fc^Wa = 0. In particular, we can 

write Ha'^Vc^^ = t.^ + , where is given by expression © [cf. Eq. (5.23) 

in Ref. [26]. Contracting with we find: l^XI ci°' = a pre-geodesic equation 

where the non-affinity coefficient is defined as n (^) = Twi^^ If a fohation {St} 
of % is given, we can write [cf. Eq. (5.35) in Ref. [26j : uj'i^ = fji^' — K^^^fca. 

Vectors l'^ and can be completed to a tetrad /c", (ei)", (62)°}, where {ciY 



®We abuse notation and employ the same notation employed in sections St of T-L, cf. Eq. l(4|l . 
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are tangent to sections St- Normalizations given in ^ are then completed to 

I ■ (e.), = , r (e,), = , (e.)"(e,)b - 5a, • (19) 

Defining the complex null vector m" = --i=[(ei)° + 1(62)"], the Weyl scalars are de- 
fined as the components of the Weyl tensor C"^^^ in the null tetrad {^°, fc", m°, m'^} 

*i = LTnH'k<^ ^-4 - mak^m'k<^ (20) 

*2 = Lm^m''k'^ 

3.1.2. Null hyper surf aces: evolution 

It is iUustrative to give a 3+1 perspective on T-L. Given a fohation % = IJteK'^* 
us evaluate explicitly the evolution along £° of quantities defined on sections St ■ 

i) Expansion equation (null Raychaudhuri equation): 

- + i^t^)' + ailV^^)"' + ^irT^.tt = . (21) 

ii) Tidal equation: 

iii) Evolution for fla'. 

S,ni'^ + = SnT^.i^q''^ + ^D, (^«(^) + ^) - '^ca(^)^ . (23) 

3.1.3. Nan- Expanding Horizons 

A NEH^Sni ig a, nuU-hypersurface H « S'^ x M, on which the expansion associated with 
vanishes (0^^^ = 0), the Einstein equations hold and —T°-c(-^ is future directed 
{null dominant energy condition). Note that any foliation % = UtGR*^* produces a 
fohation of U by MOTS 5*. 

i) NEH characterization. Making 6''^' = in the Raychaudhuri Eq. (ETI) we get 

a,b(T°^ + STrT.b^f'' = . (24) 

Since the two terms are positive-definite, they vanish independently. This provides 
an instantaneous characterization of a NEH: 

= , <j'^2 = , Tad^l"" = . (25) 

From Eq. ^ with = it follows o'-^ = 0. The NEH characterization is 
equivalent, cf. Eq. ([9]), to the evolution independence of the induced metric qab 

SiQcd - iei',^ - . (26) 

From Eq. we conclude that a NEH fixes half of the degrees of freedom in the 
second fundamental form /C^j, of St in Ai. This will be relevant in section 14.2.1 1 
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ii) Connection V a on a NEH. A null hypcrsurface has no unique (Levi-Civita) 
connection compatible with the metric. However, on a NEH T-l one can intro- 
duce a preferred connection as that one induced from the spacetime connection 
Va- u'^'^cw"' = u'^V , Vw°,it;° e TH. Indeed using NEH characterization (^5]) . 
u^Vcw" is tangent to U: Idiu"^ cw"^) = w=Vc(^dw'^) - u'^w'^Q^^} = 0. 

iii) Geometry of a NEH. We refeJ^to the pair (gab, Va) as the geometry of a NEH. 
Writing the components of the Va connection in terms of quantities on St 

kdi/cv'' = ^Daiv^k,) - q'ye['^ (27) 

the free data on a NEH are given, from an evolution perspective, by 
{qab \st , \st , K*^'' k , Q^ab \st ) i where qab is time independent. 

iv) Weyl tensor on a NEH. Under the rescaling ©, the 1-form wi*"' transforms as 
wi^^ — > Loli^ + Vain/. Its exterior derivative duj^^^ provides a gauge invariant object: 
understanding uj^f' as a gauge connection, dw*-^^ is its gauge- invariant curvature. 
Using the NEH condition, Q^^^ = 0, one can express (cf. section 7.6.2. in Ref. 1^5)) 

dw(^) = 2 Im*2^e • (28) 

Hence, Im\I'2 is gauge invariant on a NEH. Actually the full is invariant, as it 
follows from its boost transformation rules and the values of ^'o and on a NEI^^, 

*o|« = ^-ik =0 . (29) 

3.1.4. Weakly Isolated Horizons 

A Weakly Isolated Horizon (WIH) ("H, [^°]) is a NEH together with a class of null 
normals such that: SiUJa'' = 0. This condition permits to set a well-posed vari- 
ational problem for spacctimes containing stationary quasi-local horizons. This en- 
ables the development of a Hamiltonian analysis on the hor izon H leading to the 
construction of conserved quantities under WIH-symmctrieJ^. In particular, the 
expression for the angular momentum in Eq. (|15p is recovered 

with = ^Dcf (0° is an axial Killing vector, in particular divergence- free). 

The WIH structure is relevant for the discussion of IH thermodynamics (cf. A. 
Nielsen's contribution). We do not address this issue here and just comment on the 
equivalence of the WIH condition with a thcrmodynamical zeroth law. Reminding 
wi^'' = fii^'' — K^^^fca, the (vacuum) evolution equation (f23| for leads to Cifia^ = 
^ZJaK^^). More generally, S^uji^'' = Vk*^' (cf. for example Eq. (8.5) in Ref. [26]). That 
is, on WIHs the non-affinity coefficient {surface gravity) is constant: k^^-* = Kq. 
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WIHs and NEH geometry. WIHs do not constraint the underlying NEH geometry. 
In other words, every NEH admits a WIH structure. In fact, given k'^^ ^ const, the 
rescahng £' ~ ai, with Kq = const = Vfo; + aK^^\ leads to a constant k^^ ' — Kq. 
Finally, free data for a WIH are again (gabUt, ^Ist, 't'^^k, ©it^lst), but now Qablst, 
^^a^'lst and K^^'l?^ ~ Ko arc time- independent. 

3.1.5. (Strongly) Isolated Horizons 

An isolated horizon (IH) is a WIH on which the whole extrinsic geometry is time- 
invariant: [Si, Vq] = 0. This condition can be characterizecP226| g^g g^Q(k) _ fj^g^f^ 
leads to the geometric constraint 

With Eq. (pS)) . this fixes completely the second fundamental form /C^^. Free data of 
an IH, {qab\sti^a^\sti i^^^Hn = Kq), are time independent. Their geometric (gauge- 
invariant) content can be encoded in the paii0: (^i?, Im^2)- On the one hand, 
accounts for the gauge-invariant part of qab- Regarding fii^', from dcu^^^ = 2Im4'2^e 
and K^^' = const, it follows dfi'^^ — 2Im4'2%. On a sphere St we can write 
ni^^ = f^div-frco ^ ^oxact^ ^j^^^^ ^cxact ^ 2jj^g jg gauge- dependent [cf. dH])]. From 

= 2Im4'2, the gauge- invariant part of fta^ is encoded in Im4'2. 

IH multipoles of axially symmetric horizons. On an axially symmetric IH, the gauge- 
invariant part of the geometry, (^i?, Im^'2), can be described decomposed onto spher- 
ical harmonics. On an axi ally sy mmetric section St of "H, a coordinate system can 
be canonically constructec|23ESl^ such that [with A-u = 47r(i?>i)^] 

qabdx"" ® dx^ ^ {Rn f {F~hm^9de ® d0 + Fd(j) ® dcj)) . (32) 

In particular, dA = {R-u)^sm0d6d(l> {round sphere area element). We can then use 
standard spherical harmonics Yim{9), with m ~ in this axisymmetric case 

/ Yeo{e)Ye,oiO)d^A ^ {RnfSw , (33) 

JSt 

to define the IH geometric multipole^^ /„ and L„ 

In^- I ^RY„oie)d^A 

Ln = - [ lm^2 Y„o{e) d^A . (34) 

JSt 

Then, Mass Mn and Angular Momentum J„ multipoles are 

definecPEniM] 

by ad- 
equate dimensional rescalings of In and i„. 



Note the relation with the complex scalar K in Rcfs. [32] 1331 
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3.1.6. Gauge freedom on a NEH: non-uniqueness of the foliation 

Before proceeding to the dynamical case, we underline the existence of a funda- 
mental gauge freedom in the equilibrium (null) case: any foliation {St\ of a NEH 
% provides a foliation of % by marginally trapped surfaces. This is equivalent to 
the rescaling freedom of the null normal — ?> Therefore, the amount of gauge 
freedom in the equilibrium case is encoded in one arbitrary function / on St- 

Note that in this equilibrium horizon context, the relevant spacetime geometric 
object (the hypersurface H) is unique, whereas the gauge-freedom enters in its 
evolution description due to the non-uniqueness of its possible foliation by MOTS. 

3.2. Dynamical case 

3.2.1. Existence and foliation uniqueness results 

Let us introduce two fundamental results following from the application of geometric 
analysis techniques to the study of dynamical trapping horizons. 

Property 1 (Dynamical horizon foliation uniquenesl^) . Given a dynamical 
FOTHTi, the foliation by marginally trapped surfaces is unique. 

This first result identifies an important rigidity property of DHs: the uniqueness 
of its evolution description. This is in contrast with the equilibrium null case, with 
its freedom in the choice of the foliation. In particular, on a dynamical FOTH the 
evolution vector is completely determined: h°' unique up to time reparametrization. 



satisfying an appropriate stability condition on a Cauchy hypersurface S, to each 
3+1 spacetime foliation {Yjt)t£K H corresponds a unique dynamical FOTHs % con- 
taining Sq and sliced by marginally trapped surfaces {St\ such that St CI ^t- 

This second result addresses the Initial Value Problem of DHs, in particu- 
lar the existence of an evolution for a given MOTS into a dynamical FOTH. 
The resu lt requires a stability condition (namely, Sq is required to be stably 
outermos^^^^^ , so that the sign of the variation of 0*^^^ in the inward (out- 
ward) direction is under control. This is essentially the outer conditiorP^ in the 
FOTH characterization. 

3.2.2. 'Gauge' freedom: Non-Uniqueness of Dynamical Horizons 

The evolution of an AH into a DH is non-unique, as a consequence of combining 
Properties 1 and 2 above. Let us consider an initial AH Sq C Eg and two different 
3-1-1 slicings {Sj^} and {St^}, compatible with Eg. From Property 2 there exist DHs 
Hi ~ IJti ^2 = Ut2 J with Sti = "Hi n and St2 ="^2 Ef^ marginally 

trapped surfaces. Let us consider now the sections of Hi by {Et^}, i.e. St^ = "HiHEi^ , 
so that Hi = IJj^ S't^. In the generic case, slicings {S't^} and {St^} of Hi are differ- 
ent (one can consider a deformation of the slicing {Et^}, if needed). Therefore, from 
the foliation uniqueness of Property 1, sections S't^ cannot be marginally trapped 



Property 2 (Existence 




trapped surface Sq 
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Fig. 1. Illustration of the DH non-uniqueness. Dynamical horizons Hi and 'H2 represent evolu- 
tions from a given initial MOTS corresponding to different spacetime 3-|-l slicings. 

surfaces. It follows then that T-Li and T-L2 are different as hypersurfaees in A^: if 
"Hi = 7^2, sections St^ (MOTSs) and (non-MOTSs) would coincide by construc- 
tion, leading to a contradiction. In addition to this non- u niqu eness, DHs interweave 
in spacetime due to the existence of causal constrainta^^l a DH "Hi cannot lay 
completely in the causal past of another DH 'H2 (cf. Fig. [1}. 

Comparing with the discussion in section 13.1.61 on the uniqueness and gauge- 
freedom issues in the equilibrium case, we conclude from the previous geometric 
considerations that the dynamical and equilibrium cases contain the same amount 
of gauge freedom, namely a function on 5, although dressed in a different form. 
More specifically, whereas in the NEH case there is a fixed horizon, with a rescaling 
freedom — )■ / function on St), in the DH case the foliation is fixed, but a 
(gauge) freedom appears in the choice of the evolving horizon (lapse function N on 
St). In other words, in the dynamical case the choice is among distinct spacetime 
geometric objects, T-Li and H2, whereas in the equilibrium case the choice concerns 
the description (foliation) of a single spacetime geometric object T-L. 

3.2.3. FOTH characterization 

As discussed in l2.21 a FOTH with evolution vector ft," = £°- — Ck'^ is characterized by: 
i) a trapping horizon condition: 0^^^ = 0, 5h0^^^ ~ 0, ii) a future condition 0^^^ < 0, 
and iii) an outer condition; 5^0^^^ < 0. These conditions can be made more explicit 
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in terms of the variationJ22llQl 



Sal 



(I) mo-b 



(35) 



with a and j3 functions on St- Making /3 = 



1, the outer condition writes 
- i^i? + SnTabk^'i'' < 0, 



(36) 

C) is 
(37) 



- SnTabT 



for some f and whereas the trapping horizon condition (with a ~ 1, 13 
that is 

2 J 

This eUiptic condition on C, in particular through the appUcation of a maximum 
principle relying on the outer condition Jfe^^^-* < 0, is at the heart of the non-local 
behaviour of the worldtube IJteR discussed in section 12.41 

Remark on the variation/ deformation/ stability operator 5y6^^\ Before proceed- 
ing further, Eq. psp requires some explanation. In section [2.1.31 we have introduced 
Sv in terms of the Lie derivative on a tensorial object. However, the expansion 9^^^ 
is not a scalar quantity in the sense of a point-like (tensorial) field defined on the 
manifold M. The expansion is a quasi-local object whose very definition at a point 
p £ M requires the choice of a (portion of a) surface S passing through p. In 
this sense, 6jy (with 7 a function on S) cannot be in general evaluated as a Lie 
derivative. Consider a displacement of the surface St by a vector 711". The sur- 
face St+st and therefore fl^^'jt+^t depend on the angular dependence of 7, so that 
S-yvO'^^^ /= 7(51,6*'^^. The operator (5„ still satisfies a linear property for constant lin- 
ear combinations, (Jatj+b^fl^^-' = aSv9^^^ + hSwO'^^^ {a,b £ M), and the Leibnitz rule, 
<5„(76»(^)) = ((5„7)6'(^) -t- 7(5„6'(^). Details about this operator can be found in Refs. 
[551 mi H^l- Here we rather exploit a practical trick for the evaluation of (5^1,6'^^-', 
based on the remark that given the vector normal to 5, and not multiplied by a 
function on 5, it still holds formally 6y6'^^'> — Cy6'^^\ Then, we can evaluate S^yO'^^^ 
as S^yO'^^^ — SyO'^^^ — Ci9'^^\ with u° = ^v"-. In particular, the application of this 
strategy to the second line of ((35)) goes as follows. We write fc° — f3k°- and calculate 
Si^9^^^ through a Lie derivative evaluation. This results in 



(-fc%) 



SnTabky'' 



(39) 



Using (-fc%) ^ /3, ni""^ = ni^^ 

in (jHSj) follows (cf. footnote |b|. 



-2D„ln/3 and fli^^ 



-ni^^ the expression for Sj^O^^'' 



^See also the treatment in terms of Lie derivatives in the double null foliations treatment in Refs. 

MM 



August 12, 2011 0:9 WSPC/INSTRUCTION FILE JaramilloShanghai 



16 J.L. JARAMILLO 

3.2.4. Generic properties of dynamical FOTHs 

We review some generic properties of dynamical trapping horizonj- ^'^ l ^^ l ^'^ l '^^ l ^^ i 

i) Topology Law: under the dominant energy condition, sections St are topological 
spheres. This can be shown by integrating 6k9^^^ < on iSj. Under the assumed 
energy condition, the Euler characteristic x 

is positive and, being St a closed 2-surface, its spherical topology follows. 

ii) Signature law: under the null energy condition, H is completely partitioned into 
null worldtube sections (where (5^0^^^ = 0) and spacelike worldtube sections (where 
SiO'^'^^ 7^ at least on a point). Applying a maximum principle to the trapping 
horizon constraint condition, Eq. (|57)) . it follows that either C = const > 0, or C is 
a function C > everywhere on S (cf. discussion in l2.4p . 

iii) Area law: under the null energy condition, if Si6'-^'> ^ somewhere on St, the area 
grows locally everywhere on St ■ Otherwise the area in constant along the evolution. 
This follows from applying the future condition, 0'^'^' < 0, and the signature law to 

Sh^^ ^ -ce^'''> ^ [cf. Eq. (dZD]. 

iv) Preferred choice of null tetrad on a DH. According to the foliation uniqueness 
and existence results discussed in l3.2.11 there is a unique evolution vector ft," tangent 
to H and orthogonal to St, such that h"" transports St G T,t onto iSt+^f € ^t+st- that 
is, Sht ~ 1, for a given function t defining a 3+1 spacetime foliation {S*}. Denoting 
the unit timelike normal to Et by n", the lapse function by N, i.e. Ua ~ —NVat, 
and the normal to St tangent to by s", we can write on the horizon 'Hjv 

ft" = iVji" + 6s" , (40) 

for some b fixed from N and C in ([T6l) . as 2C = {b — N){b— N). The expression of the 
evolution vector as ft" £" - Cfc" [cf. Eq. (HH)] links the scaling of £" and fc" to that 
of ft". In particular, is singled out as the only null normal to St such that ft" — > ^" 
as the trapping horizon is driven to stationarity (C — > <^ 6^9^^'' — > 0). Writing 
generically the nuU normals at JIn as = f ■ (n" + s") and fc" = (n" - s")/(2/), 
Eqs. (Uni) and lead to a preferred scaling of null normals on the DH Hn 

e% = ^ in"" + s^) , k% = ^ (n" - s") . (41) 

3.2.5. Geometric balance equations 

One of the main motivations for the development of quasi-local horizon formalisms 
is the extension of the laws of black hole thermodynamics to dynamical regimes. 
This involves in particular finding balance equations to control the rate of change 
of physical quantities on the horizon, in terms of appropriate fluxes through the 
hypersurface. This is an extensive subject whose review is beyond our scope. In the 
spirit of the present discussion, we restrain ourselves to comment on the balance 



August 12, 2011 0:9 WSPC/INSTRUCTION FILE JaramilloShanghai 



LOCAL BLACK HOLE HORIZONS IN THE 3+1 APPROACH TO GENERAL RELATIVITY 17 



equations for two geometric quantities on St : the area A = dA = \ and 
the angular momentum J[(/)] in Eq. (jisp. for an axial Killing (or, more generally, 
divergence- free) vector 0" . That is, we aim at writing 



'-^^f F-dA 



dt 



= / F-' dA 

J St 



(42) 



for appropriate area F^ and angular momentum F"^ fluxes, with d/dt associated to 
the foliation Lie-transported by h"' . Eventually, one would aim at writing a 1st law 
of thermodynamics by appropriately combining the previous balance equations 



dA 
dt 



dJ[(f>] 
dt 



F^ dA , 



(43) 



for some functions and fit on St , so that F^ is interpreted as an energy 
fl^ ^41 | 24 | 42 | 43 | 44 | 22|45 | 4 5E71 ^ gj.g^ ^^^^ towards (gH) we write evolution equa- 
tions for the expansion 0^'') and the form lli^' along the evolution vector These 
equations are given by the projection of some of the components of the Einstein 
equations onto Ti. Introducing a ^-momentum current density p a = —TabT^, with 
T° the vector orthogonal to H defined in ([T5)) , such equations provide three of the 
components of Pa- The fourth is given by the trapping horizon condition (|38p . In 
brief: 



i) Evolution element of arecPllSl (p^h" = -TabT^h"): 



2'D''Qa+STTTabT''h 



a ub 



"§7 



and kC^) = -h'^k'^V 



bf-c- 



with Oa = 4^ [Cn):^ - l/^-'DaC 

ii) Evolution normal (rotation) 

fornl21149| (^^^fc^ ^ -TbcT'=q\): 
[6h + e^^^) 17^ = ^T^a^t") - 2z?Vi;) - 'DaO^"^ + 8nq\Tbcr^ O^^^^DaC . (45) 

iii) Normal component [paT^ = —TabT^T°-)'- linear combination, using = 2£"- — h°-, 
of TabT°'h^ (area element evolution) and Tabr'^i^ [trapping horizon constraint (p8|) ]. 



In order to derive the evolution equation for A, we write A ^ dA = \ so 
that, using the transport oiSt into St+st by /i", we have ^ = Jg Sh{dA) = 9^'^'>dA 
^ = Is {^hO^'''> + (6l(''))2) dA. From Eq. (04]) it then follows 

ti-KlabT n +a, a' ' H — 



and 



(PA 
dt^ 



,dA 
dt 



St 



+ (S' - k')^^''^ 



(46) 



where k' = k — InC and k' = K{t) = A ^ J^^ k'\. Note that this is a second-order 
equation for the arecPSl. Near equilibrium, the second time derivative as well as 
higher-order terms can be neglected leading to the Hawking & Hartle expressiorl^ 



,dA 

i —r- 

dt 



St 



snTabee' 



anb 



ah 



dA 



August 12, 2011 0:9 WSPC/INSTRUCTION FILE JaramilloShanghai 



18 J.L. JARAMILLO 

Regarding the evolution equation for J[4)\, we make use of Eq. together with a 
divergence-free condition on cff" (that relaxes the Killing co ndition) a nd the condi- 
tion that (jf- is Lie-dragged by the evolution vector 

with the second term in the right hand side accounting for a non-Killing cjf' . Inter- 
estingly in dynamical (spacelike) horizons "H, the conditions ^Da4>'^ = and Sh'p"' 
completely fi#Bthe form of the vector 0": (/)° = \"-'''^DbO'-''\ 

3.2.6. Open geometric issues and physical remarks 

To close this generic section on geometric aspects of dynamical horizons, we list 
some relevant open geometric problems: 

i) Canonical choice of dynamical trapping horizon. DHs are highly non-unique in a 

given black hole spacetime. A natural question concerns the possibility of making a 

canonical choice. Ther e has been some attempts in this direction based on entropic 
argument^48i53i541551?, 

A very inter estin g avenue lies on the recently introduced 
notion of the core of the trapped regio-J^^ (see also J.M.M. Scnovilla's contribution). 

ii) Asymptotics of dynamical horizons to the event horizon. One would expect DHs 
to asymptote generically to the event horizon at late times. This is indeed a topic 
of active researctPS^'^ 

iii) Black hole singularity covering hy dynamical horizons. In addition to the asymp- 
totics of DHs to the event horizon, it is also of interest to assess their behaviour 
at the hirth of the black hole singularity, in particular their capability to separate 
(dress) singularities from the rest of the spacetime (see section [5A4|) . 

DHs as physical surfaces. Dynamical horizons are objects with very interesting 
geometric properties for the study of black hole spacetimes. In addition, from a 
physical perspective it is remarkable that they admit a non-trivial thermodynamical 
description (cf. A. Nielsen's contribution). However, it is also important to underline 
that, if thought as boundaries of compact physical objects (in the sense we think, say, 
of the surface of a neutron star), then they have non-standard physical properties: 

a) They are non-unique. From an Initial Value Problem perspective, the question 
about the evolution of a given AH is not well-posed, since it depends on the 3-1-1 
slicing choice (such non- uniqueness in evolution is typical of gauge objects). 

b) Dynamical trapping horizons are superluminal, something difficult to reconcili- 
ate with the physical surface of an object. 

c) DHs show a non-local behaviour. For instance, they grow globally (reacting as 
a whole) when energy crosses them at a given local region (even a point). This 
is a consequence of their intrinsic elliptic^ rather than hyperbolic, behaviour. 
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4. Black hole spacetimes in an Initial-Boundary Value problem 
approach 

In the context of an Initial-Boundary Value Problem approach to the construction 
of spacetimes, dynamical trapping horizons play a role at two levels: i) first, as an 
a priori ingredient to be incorporated into a given PDE formulation of Einstein 
equations, and ii) as an a posteriori tool to extract information of the constructed 
spacetimes. In this section we address their application as an a priori ingredient. 

4.1. The Initial Value Problem in General Relativity: 3+1 
formalism 

Our general basic problem is the control^ of the qualitative and quantitative as- 
pects of generic solutions to Einstein equations in dynamical scenarios involving 
a black hole spacctimc. The Initial-Boundary Value Problem approach provides a 
powerful avenue to it. Such a strategy is well suited, on the one hand, to the use of 
global analysis and Partial Differential Equations (PDE) tools for controlling the 
qualitative aspects of the problem and, on the other hand, to the employment of 
numerical techniques to assess the quantitative ones. In particular, we focus here 
on the Cauchy (and hyperboloidal) Initial Value Problem. 

4.1.1. Einstein equations: Constraint and Evolution System 

General Relativity is a geometric theory in which not all the fields constitute physical 
degrees of freedom (gauge theory), so that constraints among the fields are present. 
In the passage from the geometric formulation of the theory to an analytic prob lem 
in the form of a specific PDE system, several PDE subsystems enter into scenJ^. 
First, the constraint system is determined by the (Gauss-Codazzi) conditions that 
data on a 3-dimensional Riemannian manifold must satisfy to be considered as 
initial data on a spacetime slice. The Hamiltonian and momentum constraints are 
determined by the Gabn^ components of the Einstein equation, where n" is a unit 
timelike vector normal to the initial slice. Second, the evolution system is built 
from the rest of Einstein equation, including possible auxiliary fields. The gauge 
system determines the dynamical choice of coordinates in the spacetime. Finally, a 
subsidiary system controls the internal consistency of the previous systems. 

4.1.2. 3+1 formalism 

We introduce some notation regarding the 3+1 formalisrrPH. As in section 13.2.41 
given a 34-1 slicing of spacetime by spacelikc hypcrsurfaces {St}, the unit timelike 
normal to Et is denoted by and the lapse function as N , Ua = —NVat-, with t 
the scalar function defining the 3-1-1 slicing. The 3+1 evolution vector is denoted by 
t° = Nn'^ + P'^, where (5'^ is the shift vector. The induced metric on St is denoted by 
7a&, i.e. 7a6 = gab + naUij. We choose the following sign convention for the extrinsic 
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curvature of St in A^: Kab = —j'^a^cn-b = —h^injab- In particular, we can write 



~ UikDjP + jjkDiP — jijj , where the dot denotes the derivative Ct- 
Indices i,j,k... are used for objects leaving on St. For concrcteness, we focus on a 
particular 3+1 decomposition of Einstein e quat ions, namely involving the following 
conformal decomposition {conformal Ansat^^ for data (-yij^K^^) on St 

7,, = , K,, = + ^K^,, , (48) 

for several ( choices. Denoting by Di the Lcvi-Civita connection associated with 
7ij and inserting (|48p into Einstein equations leads to a coupled elliptic-hyperbolic 
PDE system on the variables f3'\ N and jab- The elliptic part has the form 

DkD'^fi' + ^D'Dkf3'' + ^R\(3'' = Sp['f,N,p\K,^,...] (49) 
DkD^N + 2bk\n'^ b^N = S'iviiV, if, 7, A', ...] , 

where the equation on 'I' follows from the Hamiltonian constraint, the equation on 
/?' follows from the momentum constraint and the third equation on follows from 
a (gauge) condition imposed on K. If only solved on an initial slice with 7,-,-, 7'-', K 
and K as free data, this system constitutes the Extended Conformal Thin Sandwich 
approach to initial 

datJSEl. 

If we solve it during the whole evolution, together 

with 

^ - I^Af^- - 2/:^^ + CpCpf^ = 5y[iV,*,/3\if,7, ...] , (50) 
for 7y , it defines a particular constrained evolution 

formaUsnJSSESinZl. 



4.2. Initial Data: Isolated Horizon inner boundary conditions 

There are two standard approaches to ensure that initial data on a slice Eq corre- 
spond to a black hole spacctime. The punctures approach exploits the non-trivial 
topologyfi^lSS' of So, whereas the excision approach removes a sphere from the ini- 
tial slice and enforces it to be inside the black hole region. In a sense, they both 
reflect the global versus quasi-local discussion in section [1] Here we discuss the use of 
inner boundary conditions derived from the IH formalism, when constructing initial 
data of black hole instantaneously in equilibrium in an excision approach. 



4.2.1. Non-Expanding Horizon conditions 

The NEH condition gJ^^J = in Eq. ([26]) [or ([25]) ] provides three inner boundary 
conditions for the elliptic system (|49|) . In particular, they enforce the excised surface 
Sq to be a section of a quasi-local horizon instantaneously in equilibrium. 

For a given choice of free initial data in system (|49| . the geometric NEH inner 

ill) 

boundary conditions, = 0, must be complemented with two additional inner 
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boundary (gauge) conditions. Denoting by the normal vector to St tangent to Sj, 
we write /?' = P^s^ + with = /3*Sj; and /3^s,; = 0. Adapting the coordinate 
system to the horizon (i.e. = 4=> (3^ = N) supphes a fourth gauge condition 

that, together with the 6l(^' = and a^^^ = NEH conditions, re&d^^^^^^ 

^bapl+^bt,Pl-CbcPl)qab^Q , P^^N, (51) 

where qab = ^^^ah and j3\ — <7a&/3"''. A fifth boundary condition, namely for A^, can 
be obtained by choosing a slicing inner boundar y cond ition. The (gauge) weakly 
isolated horizon structure can be used in this sensi^^I^. 



4.2.2. (Full) Isolated Horizon conditions 

The next geometric quasi-equilibrium horizon structure is a (full) IH (cf. sections 
I3.1.4l and l3.1.5p . This involves three additional conditions that cannot be accommo- 
dated in system ()49p for fixed free initial data. However, we can revert the argument 
and employ IH conditions to determine improved quasi-equilibrium free initial data 
7a6 and ^ab by solving the full set of Einstein equations and ([50]) under a 
quasi-equilibrium Ansatz. Namely, we can set dtj"'^ and ^-^r- in (|50p to prescribed 
functions ff^ and and consider the elliptic system formed by together with 

- ^AT' + CpCpT' = - /2 " + 2Cpff . (52) 

This extended elliptic system is solved for ten fields: {'i>,/3°',N) and the five 7"''. 
Geometrically, we need to impose four gauge inner conditions, leaving exactl y six 
inner conditions to be fixed. Remarkably, this fits exactly the six IH condition 

©11^= , ei',)=eitH^o,<7a.,f^i'') ^ Ffr(«.,vI',r,Ar,7,,) = 0, (53) 

where -F^* ' is determined by the expression for Q^J^j^ in Eq. pip , fixed up to the 
value of the constant Kq. It is interesting to remark that this IH prescriptiorPSI 
completely fixes (up a Kq one-parameter family) the extrinsic curvature tensor /CJ^^ ~ 
k^Qab + ^"^il^ [cf- Eq. dll)] of So as embedded in the spacetime M. 



4.3. Constrained evolutions: Trapping Horizon inner boundary 
conditions 

The elliptic- hyperbolic system (|49 |) -(f50 | provides a constrained evolution scheme for 
the dynamical construction of the spacetime. Adopting a excision approach to black 
holes, we need five inner boundary conditions for the elliptic part of the system. In 
principle, dynamical trapping horizon conditions on the inner boundary worldtube 
H =^ UfSt provide a geometric prescription guaranteeing that H remains in the 
black hole region. However, imposing FOTH conditions on H can be too stringent 
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in generic evolutions. The reason is that the constructed worldtube of MOTS "H, re- 
garded as a hypersurface in spacetime, can change signature. This is in con flict with 
the outer condition in l2.2l fsomething related to jumps occurring genericalljE^I^SI^ 
in AH evolutions; see l5.1.1|) so that the resulting PDE system can become ill-posed. 
In this context, trapping horizon conditions together with the requirement of recov- 
ering NEH inner conditions at th e eq uilibrium limit, provide an appropriate relaxed 
set of inner boundary conditionJ^. More specifically, trapping horizon conditions 
provides two geometric conditions 6^^^ = and ShO''^^ ~ 0, whereas three additional 
gauge conditions guarantee the recovery of NEH at equilibrium. 

As a first step, as in 14. 2. 11 we choose a coordinate system adapted to the horizon. 
This means that spacetime evolution t° is tangent to Decomposing the shift as 

^ ^j_^o_^^a^ ^^^^^ ^j-ittcn as t° = iVn"+/3° = {Nn''+bs'')+l3^^+{l3^-b)s'' = 
h°- + l3^^ + - 6)s". Therefore t" is tangent to n if and only if /^-^ = b. 

i) Geometric trapping horizon conditions. Condition ^^^^ — leads, in terms of the 



3-1-1 quantities in 14.1.21 to the expression in the first line of Eq. ([5T|) . Condition 
J^^t^) = in Eq. dSHl), using the adapted coordinate system f3-^ = b, leads to 

[-'^Da^D'' - 2L"^Da + A] {13-^ - N) = B{I3^ + N) , (54) 

where La = K,jS''q\, A=\'^R- ^DaL" - LaL" - 47rr^^(n'' -I- s^)(?i'' - s"), and 
B = la[^^ a^^^"'' + A-KTabin'' + s^){n^ + s^), with t ^ + s". 

ii) Gauge boundary conditions I. Aiming at recovering NEH boundary conditions 
for l^ff, we first express ShQab = lab + '2.(7^^^ in adapted coordinates (/i" = — /3ff) 



2<^i6^ = - ^ In \/9 - ( '"Daf^l + ''D,Pl - ''DJil q^b ) , (55) 

Then, the coordinate choice dtQab — dt In-y^ Qab = leads to the condition on Pa 

^dJI + ^DbA - ^D,(3l qab = -2a'^al' , (56) 
that is completed by using the evolution equation for cr^^'' on H 



f dqab d 



?j - [^Da/3l + ^DbPl - ^D.Plqab) 



Ah) 
'ab 



-q\qf,C'defici' - C^q\qf,C-defkck^ 



SirC 



q\q\T,d - ^{q-''T,d)qab 



(57) 



iii) Gauge boundary conditions II. The slicing condition for N is essentially free. 
However, from Properties 1 and 2 in section 13.2.11 such a choice is equivalent to 
choosing a dynamical horizon H. Since each "H is a genuine geometric object, this 
suggests the possibility of recasting into geometric terms the gauge choice of inner 
boundary condition for N , by selecting a trapping horizon H satisfying some specific 
geometric c riterion for Ji. As an example of this, maximizing the area growth rate 
i of H leadsPIZHlto the condition p-^ - N = -const • 6l(^), with fc" ^n"- s". 
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5. A posteriori analysis of Black Hole spacetimes 

We address here the application of dynamical trapping horizons to the a posteriori 
analysis of spacetimes, their main application in the Initial Value Problem approach. 

5.1. "Tracking" the black hole region; Apparent Horizon finders 

As discussed in II. 2[ event horizons cannot be located during the spacetime evo- 
lution. However, in applications such as numerical relativity, assessing if a region 
of spacetime lays inside the black hole region can be crucial during the evolution. 
Under the assumption of cosmic censorship, the location of AHs in spatial sections 
Yit and the worldtubes constructed by piling them up (see ll.3.ip are extremely use- 
ful to determine the evolutive properties of the black hole. In this sense, apparent 
horizon finders prove to be extraordinary practical tools. These are algorithms for 
searching surfaces St C Sf that satisfy the MOTS condition 9^^^ = 0. There are 
many approaches to this problerrP^l, but all of them aim at solving the condition 
DiS^ — K+KijS^s^ = 0. For instance, assuming spherical topology, we can character- 
ize the surface in an adapted (spherical) coordinate system as F(r, 8, ip) = r — h{9, ip) 
with F = const, so that the normal vector to St is given by s, = , ^ DiF with 

^ ^D'FD.F 

DiF = {1, —deh, —diph) in the spherical coordinate system. The MOTS condition 
becomes then a non-linear elliptic equation on h that can be solved very efficiently. 

5.1.1. Understanding apparent horizon jwrnps 

Non-continuous jumps of AHs occur generically i n 3+1 bl ack hole evolutions. The 
dynamical trapping horizon framework sheds lighil^ ^'^ l ^^ l on these AH jumps, sug- 
gesting a spacetime picture where the jumps are understood as multiple spatial cuts 
of a single underlying spacetime MOTS worldtube. Jumps are associated with the 
change of metric type of the horizon hypersurface (see Fig. [2]) . This is particularly 
dramatic in binary black hole simulations, where at a given time t the two individual 
non-connected horizons jump to a common one. A specific pred iction of the dynami- 
cal horizon picture is that new (common) horizons form in pair^^^^. the outermost 
(apparent) horizon growing in area and a dual inner one whose area decrease in the 
time t. Apart from providing a better understanding of the underlying geometry 
of the trapped region, this spacetime picture can be of use in the study of flows 
interpolating between a given MOTS and the eventual event horizon, something of 
potential interest for studies of the Penrose inequality (sec below [572. 2p . 

5.2. Horizon analysis parameters 

Assigning parameters to (individual) black holes can offer crucial insight into the 
dynamical evolution. These can be physical parameters like the mass or the angular 
momentum, or diagnosis parameters informing of relevant dynamical properties. 
Given the generic absence of background rigid structures, first-principles parameters 
are often out of reach and one must follow non-rigorous or pragmatic approaches. 
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Fig. 2. Illustration of AH jumps as multiple cuts of a single spacetime MOTS-worltube H. In 
particular, timelike sections of H produce jumps (null hypersurfaces are represented with 45°). 

5.2.1. Mass and Angular Momentum. IH and DH multipoles 

In our discussion we have avoided entering into first-principles physical issues, stress- 
ing rather the geometric properties of dynamical trapping horizons and their appli- 
cations. However, mass and angular momentum estimates for individual black holes, 
either fundamental or effective, are extremely important in the modeling of astro- 
physical systems involving matter or binary systems. The problem has two aspects. 
First, one must identify a surface to be associated with the black hole boundary. 
Discussion in section [T] shows that this is a delicate question. In any case, AHs 
provide surfaces St G St tracking the black hole region, that can be employed as 
preferred choices for pragmatic estimations. The second problem refers to the ambi- 
guities in the quasi-local charact erizatio n of the gravitational field mass and angular 
momentum in General 

RelativitjlS1182| 

Regarding the angular momentum, the Ko- 
mar express ion (1151) characterizes appropriately the axisymmetric case. Effective 
prescriptionJSSlMlHS] g-jj-jgt foj- generic horizons. Regarding the mass, the irreducible 
mass il/iricd A ~ 167rM^j.j,j provides a purely geometric estimation in terms of the 
area. Its physical interpretation as the portion of the black hole mass that cannot 
be extracted by a Penrose process, together with its equivalence with the Hawking 
energy, A/nawking = v/^/(167r)(l + l/(87r) §e^'^^e^^UA) for MOTSs make s it useful 
in numerical applications and in the thermodynamical treatmentJSHlSl Given A 



and J one can also considep^^ the Christodoulou expression for the Kerr mass 

There are many prescriptions for the quasi-local masJ^HSH. It is therefore crucial to 
choose and keep consistently a prescription when comparing different solutions. In 
this latter sense, the mass and angular momentum horizon geometric multi poles / „ 



and Ln in (|34| offer a useful and refined diagnosis tool in numerical studie?^^^^ 
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5.2.2. Useful diagnosis parameters 

Insight into the geometric properties of MOTS worldtubes leads to useful diagnosis 
parameters for monitoring dynamical evolutions. Geometric black hole inequali- 
ties provide a particular avenue. In particular, the conjectured Penrose's inequality 
A < 167rAf^j^j^^ for asymptotically flat spacetimes provides a bound to the AH area 
(strictly speaking, the bound is on the area of a minimal surface enclosing the AH). 
A violation of ep„„„3„ = A/ {IGttM'^^^) < 1 indicates a more exterior MOTS. In the 
axially symmetric case this can be refined in terms of a so- 

calledSZEHl 

Dain number 

Moreover, the rigidity part of the conjecture provides an extremely simple character- 
ization of Kerr as satisfying e^^j^^ = 1. In the same spirit, the geometric inequalitjJ^ 
J < M^j^j^ provides a characterization of (sub)extremality of black holes. However, 
these inequalities involve total quantities such as the ADM mass. It is remarkable 
that the dynamical horizon structure (actually the owter trapping horizon condition ) 
provides exactly the needed conditions to prove the quasi-local inequalitj^^niSIEll 

A > 87r| J| , (60) 

in generic spacetimes with matter satisfying the dominant energy condition. The 
validity of the area-angular momentum inequality (|60p is equivalent to the non- 
negativity of the surface gravity n of isolated and dynamical horizoni^, supporting 
the internal consistency of their first law of black hole thermodynamics. Inequality 
(jnni) provides a quasi-local characterization of black hole (sub)extremality, that is 
directly related to changes in the horizon metric typj^ and jumps discus sed in 
15.1.11 This is also the context of the Booth & Fairhurst extremality paramete^^^^ 

e = 1 + — / dA SkO'^'^^ < 1 . (61) 

5.3. Heuristic and effective approaches in a posteriori spacetime 
analysis 

Hitherto we have discussed analysis tools to be applied in numerically constructed 
spacetimes, but related to sound geometric structures. However, when developing a 
qualitative understanding of the underlying dynamics, involving e.g. a comparison 
with Newtonian or Special Relativity scenarios, the available geometric notions are 
often not enough. This is manifest in astrophysical contexts requiring estimations 
for linear, orbital angular momentum or binding energies. In some cases, a choice 
must be done between saying nothing at all or rather adopting a heuristic approach. 

An example of the latter is the following heuristic proposal^ for a quasi-local 
black hole linear momentum. Given a vector transverse to a MOTS S, applying 
on S the linear momentum ADM prescription at spatial infinity leads to 



PiO = i / (Kab - K^ab) Cs' \ ■ (62) 
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In spite of its ad hoc nature, this quantity has been successfuUy applied in the 
analysiJ^ of linear and orbital angular momentum in binary black hole orbits and 
in the recoil dynamics of the black hole resulting of asymmetric binary mergers. 

5.4. An effective correlation approach to the analysis of spacetime 
dynamics 

The qualitative and quantitative understanding of strong-field spacetime dynam- 
ics represents a challenge in gravitational physics both at a fundamental level and 
in applications. In astrophysical settings a natural strategy consists in extending 
to general relativistic scenarios the Newtonian celestial mechanics approach. This 
has indeed led to fundamental achievements in the understanding of the physics 
of compact objects. However, the focus on the properties of individual objects, in 
particular in multi-component systems, also meets fundamental obstacles in a grav- 
itational theory i) without a priori rigid structures providing canonical structures, 
and ii) with global aspects playing a crucial role. The latter encompasses global 
causal issues and also the in-built elliptic character of certain objects, both as- 
pects relevant in the characterization of black holes. In this context, an approach 
to spacetime analysis that explicitly emphasizes the global/quasi-local properties of 
the relevant fields, at the price of renouncing to a detailed tracking of the geometry 
and trajectories of small compact regions, can offer complementary insights to the 
celestial mechanics approach. Such a coarse-grained effective description is much in 
the spirit of the correlation approach in the analysis of complex condensed-matter 
systems or in quantum/statistical- field theory, where the functional structure of the 
(local) dynamical fields is encoded in the associated n-point correlation functional^. 
Such an approach underlines the relational aspects of the theory, as a complemen- 
tary methodology to the isolation of the dynamical properties a compact parts of 
the system. In sum, we can paraphrase the strategy as aiming at a functional and 
coarse-grained description of the spacetime geometry, by importing functional tools 
for the analysis of condensed matter and quantum/statistical field theory systems. 

5.4.1. Cross- correlations of geometric quantities ai test screens 

The strategy outline d above is admittedly vague. We sketch now a particular 
implementatiorP^I^ of some of its aspects in a cross-correlation approach to the 
analysis of spacetime dynamics. Aiming at studying the gravitational dynamics in 
a given spacetime region TZ, we consider an outer Bo and an inner Bi hypcrsur- 
faces lying in the causal future of TZ. These hypcrsurfaces are taken as outer and 
inner boundaries of the bulk spacetime region of interest. The geometry of Bo and 
Bi is causally affected by the dynamics in TZ, so that Bo and Bi can be under- 
stood as balloon probes into the spacetime geometry. In other words. Bo and Bi 



''N-point correlation functions encode the functional structure of the local fields. A coarse-grained 
description appear as a truncation to a finite number of n-point functions. 
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r = i+ 




i 

Fig. 3. Cartor-Pcnroso diagram representing a generic (spherically symmetric) collapse and illus- 
trating the cross- correlation approach to near-horizon gravitational dynamics. 

provide test screens (they do not back-react on the bulk dynamics) on which we can 
construct geometric quantities ho and hi to be cross-correlated. Choosing causally 
disconnected screens Bo and Bi, a non-trivial correlation between ho and hi encodes 
geometric information about the common past region TZ. We can think of this as 
the reconstruction of the interaction region from the debris in a scattering experi- 
ment {inverse scatt ering p icture). Let us now restrict to the study of near- horizon 
spacetime dynamicJ^SESI jj^ g^j^ (asymptotically flat) black hole spacetime setting, 
null infinity and the (event) black hole horizon provide canonical choices 
for Bo and Bi, respectively (cf. Fig. |3]). Retarded and advanced null coordinates 
u and V provide good parameters for quantities ho and hi calculated as integrals 
on sections Su C J^"*" and Sy C £. A meaningful notion for the cross-correlation 
between ho{u) and hi{v), considered as time series, requires the introduction of a 
(gauge-dependent) mapping between u and v at J^"*" and £. We refer to this point 
as the time- stretching issue. 

5.4.2. Cross-correlations in an Initial Value Problem approach: dynamical 
horizons as canonical inner probe screens 

The adopted Initial Value Problem approach has a direct impact in the cross- 
correlation picture above. In particular, the event horizon is not available during 
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r = i+ 




i 

Fig. 4. Cartcr-Penrose diagram for the cross- correlation picture in a Cauehy IVP approach. 

the evolutiorQ. Instead, the (outermost) DH % fixed by the chosen 3+1 fohation 
stands as a natural spacetime inner boundary Bi. Ahhough any hypersurface cov- 
ering the black hole singularity could be envisaged for the present cross-correlation 
purposes, the DH % provides a natural geometric prescription. Regarding the time- 
stretching issue, the time function t defining the 3+1 spacetime slicing automatically 
implements a (gauge) mapping between retarded and advanced times u and v. Cross- 
correlations between geometric quantities at % and can then be calculated as 
standard time-series hi{t) and ho{t) (cf. Fig. Due to the gauge nature of t, the 
geometric information in quantities hi{t) and h2{t) is not encoded in their local (ar- 
bitrary) time dependence, but rather in the global structure of successive maxima 
and minima. The calculation of cross-correlations must take this into accountPSEU. 
This means, in particular, that quantities to be correlated must be scalars. 



5.4.3. Application to black hole recoil dynamics: towards DH news functions 

In the context of the study of black hole recoil dynamics after an asymmetric merger, 
let us take ho{u) as the Bondi flux of linear momentum along a (preferred) direction 



i(u)= lim — (/, (f Si) |A/'(u)P(ifi, Niu)^ ^Au')du' . m) 

du (u.r-i-oo) StT Jg^^ ^ ' 



'Regarding , a pragmatic choice in a Cauchy approach consists in substituting it by a timclikc 
worldtube of large radii spheres. However, can be kept if using a hyperboloidal foliation. 
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Here N is the news function at J^"*" , and is a given spacehke transverse direction 
to Su,r, so tliat {dP^ /du)[S] is a scalar. A natural choice^ for hi{v) would be given 
by the expression (|63)) with at substituted by some ^'o at "H. A preferred 
null tetrad on Sy is then needed, something that for DHs is provided by and fc^ 
in (|4ip . Using them in pO)). the preferred Weyl scalar '^^ is employed to construct 



1 



K''m{v)^-^j> {es,)N'^^\v) dA, with AA7(z;)- / {v')dv' . [M) 

In spite of the formal similarity between (|63| and (|64l) there is a fundamental dif- 
ference: whereas {dP^ / du)[(,] is an instantaneous flux through J^'^, this is not true 
for K^[£]{v). The function M{u) can be written in terms of geometric quantities 
on sections Su- This local-in-time behaviour is a crucial feature of any valid news 
function and it is not shared by J\fj^\v). However, it suffices to modify J\fj^\v) 
with terms completing the integrand '^q {v') to a total differential in time. Noting 
I'^al'^b^icfd^^i^ = 'i'ornamb + ^'o^^a'Tib, inspection of Eq. ((^ [actually its dynam- 
ical version with /i° instead of £"■] suggests the identification of a correct news-like 
function at H as proportional to the shear cr^^'' (see also Refs. lWllMl for the discus- 
sion of the news in quasi- local contexts). In tensorial notation, we write 

- -i £ ies.) (Aa'«A/S!'g) dA, with aC'^ = -^<J^Jl\ (65) 

where the coefficient in A/"^'^ guarantees the correct factor in the leading-term. 
This {dP^ /dv)[S,] provides a natural quantity to be correlated with {dP^ /du)[^]. 
The notation underlines the local character in time as the flux of a quantity P^ [£] , 
but no physical meaning is given to the latter. It is worthwhile, though, to remark 
the formal similarity of the monopolar part of the square of the news A/'jj^'^, i.e 

dE^ , ^ If (h\ „h , . 1 



2 



(«) = — 4> '^^'l'<h)dA = 

dv ^ ' WttJs^ ''^ IGtt 



with the expression of the fl ux of g ravitational energjEISH through a DH, in partic- 
ular with its transverse 

parlW 

The identification of cr^^^ as a news-like function 
suggests a further step, by introducing a heuristic notion of Bondi-like 4-momcntum 
flux through H. Considering the unit normal to H (f° = t"" / \/\t''ti,\ = 
(r + C/c°)/%/2C = (fen° + Ns'')/V2C), and for a generic spacetime vector r]" 



■iAn effective curvature vector' ISSl constructed from the Ric ci scalar ^R on sections <S^ of T-C provides 
an intrinsic prescription for hi{v) leading to non-trivialS3! cross-correlations with (dP^ /du)[(]. 
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has formally the expression of a Bondi-like 4-momenturr0. The flux of energy asso- 
ciated with an Eulerian observer n" would be 



dv dv IGtt y'2C 



(68) 



/2C 

dPl 

dv 



where ^ = y/l + N'^/2C. The flux of linear momentum for e TSt would be 

Near equilibrium (C 0), we have a'^Jl'a'^^^ - C on DHs [of. Eq. (H?])] so that 
expressions and are regular (0(-\/C')). Integrating ([69]) in time would lead 
to a Bondi-like counterparj] of the heuristic ADM-like linear momentum in ((62|). 

Before finishing this section, let us mention that the present discussion on hori- 
zon ne ws-like functions can be related^ to a viscous fluid analogy for quasi-local 
horizoni^IEH!_ Jn particular, geometric decay an d os cillation timescales (respec- 
tively, r and T) can be constructed on the horizorP^ from the expansion ^C^) and 
shear cr^^\ respectively related to bulk and shear viscosity terms. In the context 
of black hole recoil dynamic s, t his provides an instantaneous geometric prescrip- 
tion for a slowness paramete^^ P = T/t controlling the qualitative aspects of the 
dynamics. 



5.4.4. The role of the inner horizon in the integration of fluxes along % 

Flux integrations along % require appropriate parametrizations of "H, such as an 
advanced time v. Then, given the flux Fq{v) of a quantity Q{v), we can 

writeEl 

Qiv) = g(«o) +sign(C) / FQ{v')dv' , (70) 

Jvo 

this requiring an initial value Q{vo). However, such coordinate v is not natural in an 
Initial Value Problem approach. As discussed in lS.l.ll the 3+1 slicing {S*} induces 
a splitting of the DH into internal and external sections. The integration in (j70p 
can then be split into external and internal horizon parts (cf. Fig. [5]) 

t nt 

int , • , / /^\ I I TP \cxt (>l\j±l 



Q{t) = Qivo) + sign(C) / (Fq)™ it')dt' + sign(C) / {Fgr' {t')dt' + Res(i)(71) 
where the error Res(t) is Res(t) = sign(C) (^g)'"* {t')dt' . If the growth of Q is 

''An alternative expression would follow by using in I I67I I, instead of cr^^'crC'''''', the integrand in 

the DH energy flm l^'^ l ^^l^^l^'^ l, that would also include the longitudinal part Sla^'n'^^". 
'A related prescription for a DH linear momentum flux would be given by angular integration of 
the appropriate components in the effective gravitational-radiation energy-tensor of Ref. 1431 
™The coefficient sign{C), +1 for spacelike T-L and —1 for timelike "H, takes into account the possible 
integration of fluxes happening when timclikc sections of T-L occur; cf. Fig. [21 
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i 

Fig. 5. Illustration of the splitting of a DH into internal and external sections by a 3+1 slicing. 

understood as ultimately associated with some flow into the black hole singularity, 
the actual essential role of the horizon % would be that of capturing the associated 
fluxes. This assumes that the worldtube "H begins at the formation of the singularity. 
More complex singularity structures (as those coming from a binary merger) would 
require a more detailed analysis of this point. From this perspective, there is noth- 
ing intrinsically special about dynamical horizons: any hypersurface separating the 
black hole singularity from past null infinity (e.g. the event horizon) would be 
appropriate for fluxes evaluation. However, from a quasi-local perspective, if DHs 
are shown to cover systematically the black hole singularity (or, more generally, the 
inner Cauchy horizon), they actually provide excellent geometric prescriptions for 
such test screens (this is the motivation for the point iii) in section 13. 2. 6|) . 



5.4.5. Auxiliary test-field evolutions in curved backgrounds 

In 15.4.31 we have considered cross-correlations between different contractions of the 
Weyl tensor at distinct hypersurfaces. It is legitimate to question if such cross- 
correlations are meaningful at all, given their a priori different geometric contents. 
Let us consider the following approach to this issue: evolve, together with the gravi- 
tational degrees of freedom in Einstein equations, an auxiliary (set of) scalar field(s) 
$i without back-reaction on the geometry (i.e. test fields) and whose evolution on 



r = 
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the dynamically evolving background spacetime closely trackfj^ its relevant geomet- 
ric features. Then, the correlation approach outlined in 15.41 for a (coarse-grained) 
extraction of geometric content, can be applied directly on $i. We can paraphrase 
this approach as pouring sand on a transparent surface. On the one hand, this 
removes the ambiguity in the choice of quantities hi and ho at inner and outer hy- 
persurfaces. On the other hand, and more importantly, it also permits to extend to 
the bulk spacetime the (cross-)corrclation strategy between spacetime boundaries. 

6. General perspective 

We have presented an introduction to some aspects of quasi-local black holes in an 
Initial Value Problem approach to the spacetime construction. From a fundamental 
perspective, quasi-local black hole horizons provide crucial insights into the geome- 
try of the black hole and trapped regions and a sound avenue to black hole physics 
in generic scenarios. However, quasi-local black holes also meet challenges when con- 
sidered as physical surfaces of a compact object. We have adopted a pragmatic or 
effective approach in which quasi-local black hole horizons are understood as hyper- 
surfaces with remarkable geometric properties that provide worldtubes of canonical 
surfaces in a given 3-1-1 slicing of the spacetime. We have shown how they can be 
used as an a priori ingredient in evolution schemes to Einstein equations, where 
they provide inner boundary conditions for black hole spacetimes. Then we have 
illustrated their use as a posteriori analysis tools tracking and characterizing quasi- 
locally the black hole properties and providing, through their rigidity properties, 
excellent test-screen probes into the near-horizon black hole spacetime geometry. 
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